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Abstract. We present a perturbative approach for studying inflation models 
with soft departures from scale free spectra of the power law model. In the 
perturbed power law (PPL) approach one obtains at the leading order both 
the scalar and tensor power spectra with the running of their spectral indices, 
in contrast to the widely used slow roll expansion. The PPL spectrum is 
confronted data and we show that the PPL parameters are well estimated 
from WMAP-7 data. 
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1. Introduction 

Observations of the anisotropy and polarization of the cosmic microwave background 
(CMB) determine the parameters associated with cosmological structure formation to 
unprecedented accuracy [HEilS]. Inflation remains the best motivated and predictive 
early universe scenario that is invoked to specify the spectrum of initial perturbations 
for structure formation [4J. The type of perturbation that is required results from 
the simplest class of models that predict Gaussian, adiabatic, nearly scale-invariant 
perturbations Generally a single component inflaton field is considered which 

slowly rolls along a sufficiently flat potential leading to the exponential expansion of the 
universe. A simple and completely analytically tractable regime of inflation corresponds 
to the uniform acceleration approximation [6] that correspond to power law inflation 
model. The initial scalar and tensor perturbation spectra are scale free (power law) and 
are parametrized by two numbers : the scalar perturbation amplitude and a common 
spectral index. 

For the power law inflationary models [7J [8], the scale factor a(t) = t p where 
p > 1. The slow- roll parameters for this power law case are e = —5 = -. Also it is 

characterized by the (i) uniformly accelerated expansion for which, —-Ms — const and 
(ii) the perturbations in the inflaton field 5<p are equivalent to a massless scalar field 
perturbations with the effective mass squared, m 2 e ff = 0. The power law spectrum is 
given by V(k) = Ak^ 3 ' 1 * 1 and for this scale free case the scalar spectral index n s = 1. 
Here scalar and tensor perturbation spectra are of identical shape with constant spectral 
indices n s and n t = n s — 1; and the ratio of amplitudes of the tensor to scalar power 
spectra, r = — 6.2 ut [9] 

It is well known that, although primordial fluctuations spectra expected from 
inflation are likely to be approximately flat or scale-invariant (n s (k) = 1 + 
d\nV(k)/d\nk ~ 1), exact scale invariance (n s = 1) is achievable only for a very specific 
class of models [TO] , while a slightly red spectrum {n s < 1) with small running appears 
to be a generic prediction of the simplest viable one-parameter family of inflationary 
models. 

In this article, we describe a perturbative formalism for studying inflation models 
with soft departures from scale free spectra. The details of our methodology to obtain 
the primordial power spectra is discussed in Section El And the parameter estimation 
with the perturbed power law inflationary spectrum is presented in Section [31 

2. Soft deviations from scale free spectra - Perturbed Power Law Model 

Based on the ' Hamilton- Jacobi formulation' [TT1 H2] . we parametrize the inflationary 
phase by the Hubble parameter H(<fr), expressed as a function of the inflaton field </>. 
For a single inflaton field, the determination of H{4>) can be directly translated to an 
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estimation of the inflaton potential, 



V(cf>) 



8tt yv> 



1 - 
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dlnH 



(1) 



using the 'reduced Hamilton- Jacobi equation' [61 [13] . 

The field equation for the modes of the inflaton perturbation, 5<pk are given by, 
601 + 2-5<j>' k + (k 2 + a 2 m 2 eff )5<p k = , 



(2) 



where, a 2 m 2 ^ = \ — 4- and z = a 2 ^j. The effective mass can be expressed in terms of 
slow-roll parameters as p2 



m 



eff 



H 2 



and the slow-roll parameters e and 5 are given by, 
e = — — I — r | and 



(3) 
(4) 

1 for notational 



where the m-p is the Planck mass. In the rest of the paper we set 
simplicity. 

Starting from the equation of motion of the inflaton field and using the Hamilton- 

9l »J I ) ^ s a constant for the power law models 



Jacobi formalism we can show that -4- 



H 

with uniform acceleration. Thus it is clear from the Eqs. ([3]) and (j4j) that —^r- will 
be zero for the power law models. For the power law models, the scalar and tensor 
perturbation spectra are parametrized by a common spectral index v = n s — 1 = n t 
given by, 

2 

o / dlnH ' 

-2e 



v 



dlnH 

d<j> 



(5) 



We are interested in considering the small deviations from uniform acceleration in terms 
of a small d 2 \nH / dcj) 2 perturbation [15]. In the so called Perturbed Power Law (PPL) 
approach, the predicted scalar and tensor spectra are perturbed from the scale free form 
but, at the leading order maintain a constant difference between their spectral indices. 
We consider that the expansion is locally modeled to be power law with p varying 
slowly with time. Also in terms of the conformal time n we can write the scale factor 
as a{rf) = ^( 1 / 2 ~^ so that ji — 3/2 = l/(p — 1). Thus we have p, = § + and also 

we get, 

(p - 1/2) 2 z" 



H\<P) 



rj 2 a 2 (n) 



1 

r/ 2 



/4 



(6) 



where p 2 s = y p 2 — (p — 1/2) 



2 m iff 
H 2 ■ 
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Here 5<p is equivalent to the scalar field perturbations with ~ const (weakly 
depend on rj). In case of the PPL models, the spectral indices, n s ^ n t + 1. If we denote 
fi, - 1 = !/, we can derive, 



,/lu l(\ 2 _|_ / i> 2 In // 



l-e{rj) 1 — e(?7) 1 ^ _ f dh>H 



(7) 



where, 



X = ^ = -(W)- (8) 



'<9 2 lniT 

Denoting n t as v t and defining the difference between the spectral indices as, v st = v s — v t 
we get, 

' d 2 lntf" 



considering terms up to the first order corrections in perturbation. It is possible to solve 
for H((f)) in an exact form and consequently its (^-derivatives in terms of v s and v st . 
By definition and using Eq. the scalar perturbation spectrum is obtained as, 

AjfMr^s) (*) 2 J_ (10) 
sV ; 2vr \2tiJ e(rj) v ; 

and similarly the tensor perturbation spectrum, 

where A(x, y) = ^yjfv- 1 ' ^(y) is the Euler Gamma Function, ht and /ig are defined 
below the Eq. (jSJ) with the identification \i T — > ft. 

The standard slow-roll approximation does not consider the higher order 
perturbative correction terms and thus one has \it = Us = A 4 an d there the above 
expressions reduces to those in [16]. It is also interesting to note that our PPL power 
spectra expressions are similar to those obtained in [17J by considering the first order 
corrections in slow-roll parameters. With the PPL approach one can study both the 
scalar and tensor power spectra and the running of their spectral indices also, where 
as in the standard procedure one has to consider each of them independently. While 
in principle we can go to higher order PPL corrections, but we restrict up to the first 
order PPL deviation from scale free spectra that the currently available data can reliably 
capture. 

For a given power law model, the scalar and tensor power spectra can be calculated 
theoretically once we know u s and u st . Then it is interesting to determine how well 
one can constrain H(4>) using CMB observation data. The data from WMAP can be 
used to measure the leading order deviations from power law spectra, quite accurately. 
Parameter estimation with the perturbed power law inflationary model is discussed in 
the next section. 
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3. Parameter estimation with the perturbed power law model 

For Power law model, given v s and v st the scalar and tensor power spectra can 
be calculated theoretically. Then one has to perform a Markov Chain Monte-Carlo 
sampling of the parameter space to estimate the constraints on the inflationary 
parameters and the various background cosmological parameters. We make use of 
the publicly available CosmoMC package [18J, which in turn uses the CMB anisotropy 
code CAMB [19] to generate the theoretical CMB angular power spectra, C;s from 
the primordial scalar and tensor power spectra. We make appropriate modification to 
CAMB power spectrum module to incorporate PPL power spectra. For our analysis, 
we confront the theoretical C/S with the WMAP seven year data set. We make use of 
the publicly available WMAP likelihood code from the LAMBDA web site [2D] . 

Flat ACDM is considered as the background cosmological model. The priors that 
we set for the background cosmological model are listed in Table 1. f^/i 2 is the physical 
baryon density, Q c h 2 is the physical cold dark matter density, 6 gives the ratio of the sound 
horizon to the angular diameter distance at decoupling and r is the reionization optical depth. 
The pivot point is set at fco = O.OSMpc . Table 1 also gives the priors set for the inflationary 
parameters of the PPL case, where we denote n s — 1 = v s and v s — vt = v s t and A s is the 
amplitude parameter. 



Parameter 


Lower Limit 


Upper Limit 


Q b h 2 


0.005 


0.1 


Q c h 2 


0.01 


0.99 


9 


0.5 


10 


T 


0.01 


0.8 


Vs 


-0.15 


0.0 


V s t 


-0.06 


0.06 


log[10 10 A s ] 


2.0 


4.2 



Table 1. Priors for the background cosmological parameters and the PPL inflationary 
parameters 

For multiple chains the CosmoMC code computes the Gelman and Rubin (variance of 
chain means)/(mean of chain variances) R statistic for each parameter. The program also 
writes out the value for the worst eigenvalue of the covariance of the means, which should be a 
worst case. This R — 1 statistic is also used for the stopping criterion when generating chains 
with MPI. We set MPI Converge Stop parameter, that can be used to stop the chains at 0.03 
for our runs and got the worst eigenvalue: R — 1 = 0.0554. 

The best fit values, the mean (of the posterior distribution of each parameter) and 1-cr 
deviations for the background parameters as well as the PPL parameters obtained are listed in 
Table 2. And for a comparative study, the similar list for the power law (PL), power law with 
tensor (PLT), power law with running (PLR), power law with running and tensor (PLRT) 
models are also given in Table 3. Obviously, the best fit values that we obtained matches with 
those quoted by the WMAP teamfj] [21] ■ It is clear from the tables that the parameters of 

f WMAP team set the pivot point at k = 0.002 Mpc^ 1 . 
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PPL Parameter 


Best fit 


Mean & \-a 


Q b h 2 


0.0222 


n noo^+o- 0005 
u.uzzo_ . 0006 


Q c h 2 


0.1101 


1 1 n7 +0 - 0053 

u - iiu ' -0.0053 


e 


1.038 


i n39+ 003 

i.uoy_ 002 


T 


0.0881 


0.08761°,;°™ 


v s 


-0.0367 


-o.o4oo±g:8iSg 




-0.0166 


0.0027±g;8i« 


log[10 10 A s ] 


3.103 


9 an^ - 460 

z,.yuo_ 368 


A S z 


1.449 


1 024+ - 975 




0.7340 


0.7320lffi 


Age/Gyr 


13.808 


13.784±g:ig 


Vl m 


0.2660 


0.2680l°;°274 


Z re 


10.595 


10.507+1 -^1 


Ho 


70.551 


70.663 + |^ 



Table 2. The best ht values, the mean and 1-a deviations for the various input and 
derived parameters obtained by comparing PPL spectrum with WMAP7 data 



Model 


PL 


PLT 


PLR 


PLRT 


Parameter 


Best fit 


Mean & 1-a 


Best fit 


Mean & 1-a 


Best fit 


Mean & 1-a 


Best fit 


Mean & 1-a 


n b h 2 


0.0221 


u,u " J -[).()0()5 


0.0222 


u.uzzo_ () 0005 


0.0213 


0216 +U ' ,),),)8 


0.0214 


0216 +U ' ,),),te 


n c h 2 


0.1090 


1106 +(m)M 

U.11UD_ () 0Q54 


0.1116 


1106 +,UU)b3 

U.11UD_ () 0053 


0.1202 


0.1180_ 0076 


0.1175 


0.1184_ ()0g5 


6 


1.038 




1.038 


1.039™ 


1.037 


l.UOO_Q QQ2 


1.038 


1.038_ 003 


T 


0.0840 


0.0867±g;^ 


0.0829 


0.0864t™ 


0.0900 


0914 +U - UUVV 
u.uyi4_ ()067 


0.0961 


O.U914_ 006g 


n s 


0.9556 


oW™ 


0.9608 


0.9620+8;^ 


0.9034 


u - yi ' °-0.0358 


0.9108 


9146+ iyt)m 

u.yi4D_ 03g6 


n t 


NA 




0.3881 


-0.00821™ 


NA 




-0.2321 


0.0047l , ( ;;^^ 


nrun 


NA 




NA 




-0.0387 


-0.03281™ 


-0.0367 


-0 0351 +U ' U2ti5 
u.uaoi_ 0265 


r 


NA 




0.0175 


1134 4-u.na4 


NA 




0.0177 


u.ijuu_ 0343 


log[10 w A s ] 


3.050 


°- u ' 1 -0.055 


3.069 


3 072 +U.U37 
J - u ' z -0.033 


3.098 


o.uyi _o.040 


3.091 


3 099+ u (J4i! 
•5.uyy_ ()42 


A SZ 


1.935 


1 058 +U ' 94y 

1 - uo °- 1.058 


1.025 


1.0681^ 


1.411 


1-059+iX 


1.352 


1.035+-°° 


n A 


0.7381 


n „„„_+(). (JLJV8 
u - ' dzo -0.0272 


0.7270 


0-73231™ 


0.6724 


0.68571™ 


0.6910 


0.68251--- 


Age/Gyr 


13.823 


13.785+;;-.!^ 


13.814 


i3.788+;; : ;- 


13.984 


13.92lj#}$ 


13.939 


13 q9q+" 1 t>t> 
ro.yzy_ 163 




0.2619 


2675 +(Mmb 
U - ZD '°-0.0744 


0.2730 


9677+ lUm3 
u.zot i _ 0284 


0.3276 


U.J14J_ ()473 


0.3090 


U.,51(0_q 052 g 


Z r e 


10.248 


10.428ti'^ 


10.197 


10.408+i'}^ 


11.317 


11.2391-- 


11.266 


11.286L-f 9 | 36 


Ho 


70.775 


70-694^:406 


70.020 


70.663l-|- 


65.740 


67.083l-»» 


67.043 


66.903tg;J*i 



Table 3. The best fit values, the mean and 1-a deviations for the various input and 
derived parameters obtained by comparing power law (PL), power law with tensor 
(PLT), power law with running (PLR), power law with running and tensor (PLRT) 
models spectra with WMAP7 data 

The least square likelihood parameter Xeff f° r ^ ne perturbed power law (PPL), power law 
(PL), power law with tensor (PLT), power law with tensor obeying the consistency relation 
n t = — r/6.2 (PLTC), power law with running of scalar spectral index (PLR), power law 
with running and tensor (PLRT) and the power law with running and tensor obeying the 
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consistency relation (PLRTC) models are given in Table 4. The first raw gives the xlff va l ue 
obtained by fitting the power spectrum obtained for a power law model with the perturbation 
approach (PPL), to the WMAP7 spectrum. We can see that the values of first four runs are 
almost the same (the best one being that of PPL), where as we get a slightly better xlff f° r 
a model with the running of the spectral index that involve more parameters. The second 
column of the table gives the number of extra inflationary parameters compared to the simple 
power law model. For PL the inflationary paramters are the scalar spectral index n s and the 
scalar power amplitude A s . To obtain the tensor power spectrum (PLT) we need two more 
paramters, the tensor spectral index nt and the ratio of scalar and tensor power amplitudes 
r = At/As- For PLTC model, we impose the consistency relation n< = — r/6.2 also. For 
PLR model the running of scalar spectral index n run is also considered. Similarly, the extra 
paramters of the PLRT model are nt, r and n run . And r and n run for the PLRTC model 
where the consitency relation is also imposed. It is interesting to note that in the case of PPL 
just one extra parameter v st capture the role of n t , n run and in addition the running of tensor 
spectral index also, which is not possible to be estimated by the standard PL spectra as it is 
too small. 



Model 


No. of extra 
Parameters 




PPL 


1 


7474.96 


PL 





7475.14 


PLT 


2 


7475.00 


PLTC 


1 


7475.07 


PLR 


1 


7473.74 


PLRT 


3 


7473.58 


PLRTC 


2 


7473.61 



Table 4. The best fit likelihood xtff values for the perturbed power law (PPL), 
power law (PL), power law with tensor (PLT), power law with tensor obeying the 
consistency relation nt = —r/6.2 (PLTC), power law with running of scalar spectral 
index (PLR), power law with running and tensor (PLRT), power law with running and 
tensor obeying the consistency relation (PLRTC) models. 

ID marginalized posterior distribution of the background cosmological parameters 
obtained by PPL fitting to WMAP7 data is plotted in Figure 1. It follows that there are 
no significant changes in derived values of the cosmological parameters in comparison with the 
results obtained by the WMAP team assuming a power-law model of the primordial spectrum. 

Figure 2 gives the ID marginalized posterior distribution of the PPL inflationary 
parameters. The marginalized probability for v s peaks at —0.04 (corresponding to n s = 0.96 
at the pivot point), and the marginalized probability for v st is found to be maximum at 0.0027. 
The 95% marginalised limit of v s t is found to be in the range —0.016 < v s t < 0.022. The plot 
of joint 2D marginalized posterior distribution of v s and v s t is given by Figure 3. The colour 
of the figure shows how many times the CosmoMC has probed a specific area of the parameter 
space and the inner and outer closed contour lines indicate the 1 and 2-a likelihood contours. 

The best fit Cgs for power law with running and tensor (PLRT) and power law with 
perturbative method (PPL) are presented in Figure 4. The WMAP7 binned data with related 
error bars are also plotted for comparison. It is clear that the PPL spectrum gives a very good 
fit to the observed data. 
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Figure 1. ID marginalized posterior distribution of the background 
cosmological parameters (both input and derived) obtained by PPL fitting to 
WMAP7 data. Y-axis gives the probability from to 1. 



4. Conclusions and discussion 

The high quality CMB data that has become available over the past few years and other 
observational data are just approaching the level of accuracy necessary to detect deviations 
from exact scale invariance and to distinguish between different inflationary models. The data 
indicates that the departure of the spectral index from exact scale invariance is likely to be 
small, \n s (k) — 1| <C 1, which is in good agreement with predictions of the simplest inflationary 
scenarios. A perturbative procedure for studying inflation models with soft departures from 
scale free spectra is discussed in this paper. In [J5] one of the authors forecast how well one 
can constrain H(cf)) using WMAP and Planck data. The expected l-cr error ellipses in the 
(3 In H/dcf)) 2 - (d 2 In H/dcf 2 ) plane for a fixed target power law model, are now obtained with 
actual WMAP data as given in Figure 3. 

The perturbed power law spectrum is confronted with the 7-year WMAP data. The best 
fit values that we obtained for both the input and derived parameters matches very well with 
those quoted by the WMAP team for power law inflationary models. Also it is interesting to 
note that the x^s va l ue we obtained by fitting the PPL spectrum to the WMAP7 spectrum 
is very slightly better than that for a simple power law spectrum (with and without tensors). 
We find that the parameters of inflation can be well determined by this perturbative method. 
In case of the standard power law spectrum one has to consider the scalar and tensor spectral 
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Figure 2. ID marginalized posterior distribution of the PPL inflationary 
parameters v s and v s % obtained by fitting to WMAP7 data. The 95% 
marginalised limit on the parameters are shown by the vertical lines. 
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Figure 3. The joint 2D marginalized posterior distribution of v s and v s t 
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Figure 4. Best fit Cgs for power-law (PLRT) (blue, dotted), power law with 
perturbative method (PPL) (red, solid) and the WMAP7 binned data with 
related error bars 



indices n s , n t and the running of scalar spectral index n run as separate independent parameters. 
Here, in the case of PPL just the two parameters v s and v s t will capture all these. In addition, 
PPL spectra takes care of the running of tensor spectral index also, which is not possible 
to be estimated by the standard PL spectra as it is too small. In conclusion, the perturbed 
power law model helps to determine how the general models of inflation are best studied as 
'departures from the power law model'. The probable values of the 'deviation parameter' v s % 
is found to be consistent with zero (contained between —0.04 < v s t < 0.04). More precise 
data from Planck [22J would be able to measure the deviations from power law spectra quite 
accurately and our PPL approach would help to analyse and understand the results in a better 
way. 
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